Induced P-wave superfluidity within full energy-momentum dependent Eliashberg 
approximation in asymmetric dilute Fermi gases 
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We consider a very asymmetric system of Fermions with an interaction characterized by a positive 
scattering length only. The minority atoms pair and form a BEC of dimers, while the surplus 
fermions interact only indirectly through the exchange of Bogoliubov sound modes. This interaction 
has a finite range, the retardation effects are significant and the surplus fermions will form a P-wave 
superfluid. We compute the P-wave pairing gap in the BCS and Eliashberg with only energy 
dependence approximations, and demonstrate their inadequacy in comparison with a full treatment 
of the momentum and energy dependence of the induced interaction. The pairing gap computed 
with a full momentum and energy dependence is significantly larger in magnitude and that makes it 
more likely that this new exotic paired phase could be put in evidence in atomic trap experiments. 
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I. Introduction 

It is still largely an unsettled issue what happens to su- 
perfluidity in asymmetric fermionic systems. The asym- 
metry is measured by the ratio of the minority versus 
majority number densities x = nijn^. If the only inter- 
action between the fermions is between the two species 
and if that interaction is attractive, the two extremum 
cases are clear: a fully paired phase for x — \, spanning 
the full range from a Bardeen-Cooper-Schrieffer (BCS) 
superfluid to a Bose-Einstein Condensate (BEC), and 
a non- interacting fully polarized Fermi gas for x = 0. 
What happens for intermediate values of the asymmetry 
X is still a matter of debate. In the case of a unitary 
Fermi gas strong theoretical arguments [H point to the 
existence of an exotic new phase, a supersolid, with a 
airing mechanism conjectured more than forty years ago 
for asymmetries 0.25 < x < 0.75. For larger asym- 
metries (smaller x) many authors have argued that the 
system is a normal Fermi system Q , a fact which has not 
been reconciled yet with the Kohn-Luttinger conclusions 
The theoretical arguments are based largely on bi- 
ased Monte-Carlo calculations, which by design were not 
able to reveal the existence of other phases. At the same 
time the experiments Q were either insensitive or not 
designed to reveal the presence of other phases. There- 
fore the question of superfluidity of unitary asymmetric 
Fermi systems at very large asymmetries (small x) is still 
an open question, both theoretically and experimentally. 

Theory shows that away from unitarity for large asym- 
metries and zero temperature a fermion system is super- 
fluid 3. In particular, on the BEC side of the Fesh- 
bach resonance, the ground state of a system is a mix- 
ture of two coexisting superfluids: a bosonic superfluid 
made of s-wave dimers and a P-wave fermionic superfluid 
made out of the surplus majority fermions. These surplus 
fermions, which all have the same spin, interact with each 
other only through the density fluctuations of the bosonic 



superfluid in the p-wave or higher odd orbital states. P- 
wave fermionic superfluidity has been put in evidence so 
far in one system only, liquid '^He 0], but it is expected to 
appear in many other systems as well: among neutrons 
in neutron stars at densities slightly above nuclear densi- 
ties [8| , in asymmetric cold gases with S-wave zero-range 
interaction [^, Q , and in electronic systems [13] ■ There 
are many studies suggesting the creation of P-wave su- 
perfluidity in cold fermionic gases 1J| using a P-wave 
Feshbach resonance [12| (instead of a S-wave Feshbach 
resonance) , which is a different mechanism from the one 
considered by us here. The shorter lifetime of P-wave 
Feshbach molecules could prove to be a difficult exper- 
imental challenge to overcome in order to create these 
types of superfluid systems. On the other hand the S- 
wave Feshbach molecules are rather long lived and are 
routinely created in many laboratories [l3|. The argu- 
ments for the existence of P-wave superfluids in our sys- 
tem based on the BCS formalism in the weak coupling 
limit was presented in Refs. [^[^ : When all the minority 
atoms bind into dimers (for positive scattering lengths) 
and form a BEC condensate at zero temperature, the 
surplus majority atoms will interact with each other only 
through the exchange of Bogoliubov sound modes of the 
BEC dimer condensate and their pairing is allowed only 
in a P-wave or other odd partial waves, due to the exclu- 
sion principle. In this mechanism, momentum and en- 
ergy exchanged between the fermions can be comparable 
or large when compared with the Fermi momentum and 
energy of the surplus fermions. For this reason we ex- 
pect that neither the BCS approximation, nor the more 
involved Eliashberg approximation with only energy de- 
pendence 1J| , in which only the retardation effects were 
accounted for, will be sufflcient to describe quantitatively 
this kind of P-wave pairing mechanism, even though the 
calculation based on the BCS approximation in Ref [i^ is 
expected to be qualitatively correct. 

This paper is organized as follows : In section II, we 
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describe the theoretical framework and give a brief de- 
scription of our numerical calculation. In Section III, our 
numerical results on the induced P-wave pairing gap in 
homogeneous systems are discussed . In section IV, we 
discuss the induced P-wave pairing correlations in atomic 
traps and summarize our findings. 



II. Induced interaction and self-consistent equations 

The mechanism for P-wave superfluidity in asymmet- 
ric Fermi gases is very similar to the BCS mechanism of 
superconductivity for electrons. In an asymmetric Fermi 
gas on the BEC side of the Feshbach resonance (where 
the scattering length a > 0) all the minority fermions 
form Nt = Ni bosonic dimers with a ground state energy 
£b — —h'^/ma'^. The surplus Nf = — Ni fermions (all 
with spin up in our convention) have the same spin and 
they interact with each other only via the excitation of 
the Bogoliubov sound modes in the BEC condensate of 
dimers, with which they coexist 6]. The form of this in- 
teraction has been derived a long time ago by Bardeen, 
Baym and Pines [isj and in the momentum-energy rep- 
resentation has the form 



Uind{q,(^) = 



- EqiSq + 2niyUbb) ' 



(1) 



where q and uj are the momentum and energy exchanged 
between the two interacting fermions. Ubb and Ubf are 
the boson-boson and the boson-fermion couplings, re- 
spectively, nb — NJV is the boson (dimer) number 
density and Sq — q^l{2mb) is the kinetic energy of a 
dimer with mass vrib = 2m. Ubb ~ 1.27ra/m and Ufb ~ 
3.547ra/m are obtained from the fermion-fermion scat- 
tering length a > using few-body theory techniques, 
see Ref. [6[ for details. In the usual BCS approximation 
one neglects the energy dependence of this induced inter- 
action and then Uind becomes a Yukawa type exchange 
potential with a radius given by the coherence length of 
the BEC of dimers. The form of this induced interac- 
tion is generic for the coupling of fermions with sound 
modes. Since the fermions couple to density variations 
each vertex is proportional to a gradient (and thus con- 
tributes q in momentum representation) and a coupling 
constant, and the propagating sound wave is described 
by the denominator in Eq. ([T]) . While we expect that for 
large values of the scattering length the coupling strength 
Ufb could change from its weak coupling limit, the de- 
pendence of induced interaction Uind{q,^) on the energy 
and momentum should remain qualitatively unchanged. 

It was Eliashberg lj| who showed that retardation ef- 
fects are very important. Depending on the magnitude 
of Lo the induced interaction ([1]) can change its charac- 
ter from attraction to repulsion. Tolmachev has shown 
a long time ago that the role of repulsion is somewhat 
diminished in the gap equation, see Ref. [l6| . In a more 



modern language that is translated into the fact that the 
role of an attractive interaction is enhanced and a re- 
pulsive interaction is suppressed by the renormalization 
group flow. The derivation of the Eliashberg gap equa- 



tions [ij] relied on the famous Migdal's theorem [17|, 
that the vertex corrections are small in the case of elec- 
trons. The role of the vertex corrections has been ana- 
lyzed later and found to be important in the case of high 
Tc superconductors and one can expect that these 
corrections are also important in the case of dilute gases. 

In the simple Eliashberg approximation the pairing 
field has only an energy dependence. In dilute Fermi 
gases however, as we will demonstrate here, both energy 
and momentum dependence of the induced interaction 
([T|) are crucial, as we will show by comparing the solu- 
tions of the simple Eliashberg equations and of the Dyson 
equations for the fermion self-energy with full momentum 
and energy dependence. The full energy-momentum de- 
pendence of the pairing gap equations has been inves- 
tigated in condensed matter, nuclear physics and color 
superconductivity, see Refs. UJ. 

As our results show, the pairing gap values obtained 
by taking into account the full energy-momentum depen- 
dence are enhanced considerably, especially towards the 
Feshbach resonance. Thus the prospects of putting in 
evidence for the first time this new type of superfluid- 
ity, when a bosonic superfluid coexists and leads to the 
formation of a fcrmionic P-wave superfluid, are greatly 
increased. 

Here we solve the Dyson equation for the full propa- 
gator G{p) 



(2) 



where p = (lo,p) is the fermion 4-energy- momentum vec- 
tor, e{p) = p^/2m — n, n the chemical potential, rj an 
inflnitesimal small positive constant, and T,{p) is the self- 
energy. In this case the self-energy has the form 



I](p) - [l-Z{p)]uj + x{p)r3 + cj,{p)T,, 



(3) 



where Z{p), x(p)j and (p^p) are the wavefunction renor- 
malization, the Hartree-Fock potential, and the pairing 
field, respectively, and ri_3 are Pauli matrices. From 
these equations one obtains the propagator G{p) 



G{p) = 



Z{p)uj + £(p)t3 + (f){p)Tl 

[Z{p)ujY - E^{p)+ir) 



(4) 



where e{p) — e[p) -t- x{p) and E{p) — \Je^(j)) -I- l0(p)|^. 
The self-consistency is achieved through the equation for 
the self-energy 



-j7^V^nd{p - q)TzG{q)Tz. 



(5) 



We will solve these equations by performing first a Wick 
rotation w — > zw, which will make all the integrands well 
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behaved. One can show that this Wick rotation is equiv- 
alent to taking the zero temperature hmit of the finite 
temperature Matsubara equations. 

Below we will compare the solutions for the pairing gap 
by solving both the simple Eliashberg equations and the 
Eqs. ([3]l5|) for the self-energy retaining the full energy- 
momentum dependence. We will make ansdtze Z{p) oc 
^oo(p), x{p) « ym{p) and 4>{p) cx Yim{p), where Yi„,{p) 
are spherical harmonics of the unit vector p ~ pI\p\- 
Also we will limit the magnitude of the exchanged ener- 
gies |a;| < wd = chko and momenta \p\ < hku, where the 
Debye wave vector is defined by the average dimer-dimer 
(inverse) separation ko = t^ti^J^ and = nhUbb/rrih is 
the speed of the Bogoliubov sound modes. The phys- 
ical pairing gap is given by A(p) — (f>{p)/Z{p), eval- 
uated at a; = and |p| = hkf. The pairing gap 
A(p) — AimYi„i{p) will be given in units of the Fermi 
energy of the surplus fermions £/ = fi,^fcy/(2m), where 
n/ = (iV| - Ni)/V = fc3/(67r2). The effect of the wave- 
function renormalization is a change in the density of 
states and a reduction of the pairing gap. The fact that 
the equation for the self-energy is of convolution type al- 
lows us to use Fast Fourier Transform (FFT) techniques, 
after performing the integration over the angles analyti- 
cally, and solve the self-consistent equations by iterations. 
When determining the self-consistent solution of Eqs. ^ 
[5|), the value of the chemical potential was fixed so as to 
have a prescribed ratio rib/nf. 

The properties of this system are controlled by two 
dimensionless parameters, one the dimensionless strength 
of the induced interaction 



mkfribUj^ Ub 2 
— ^ — (fc/a)^ 



(6) 



and the other the ratio of the speed of sound to the Fermi 
velocity 



2nbUbb rib., . 

s = 2 ~ — 



(7) 
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FIG. 1: (Color online) In these panels the ratio of dimer versus 
the surplus fermions varies as j/ = n(,/n/ = 2, 1, and 0.5 from 
the upper panel down. The solid black line is the pairing gap 
computed using the BCS approximation as in Ref. 3. The 
solid (red) and dotted (red) line with circles are the pairing 
gaps computed with full momentum and energy dependence, 
while the solid (blue) and dotted (blue) line with squares are 
pairing gaps obtained in the simple Eliashberg approximation. 
The solid line is for the pairing gap Aio with (l,m) — (1,0), 
while the dashed line is for An with {l,m) — (1,1). Note 
that in the Eliashberg approximation An > Aio, while in the 
case of the BCS approximation there is no dependence of the 
pairing gap on the magnetic quantum number. The pairing 
gap is evaluated at zero energy u — and on the Fermi surface 
\p\ = hkf . In all cases we plotted the magnitude of the pairing 
gap, after the angular dependence Yim{p) has been factored 
out. The region of smaller values of fc/a is shown in inserts. 



The asymmetry parameters x and y = nb/nj are related 
y = x/{l — x). In defining the strength of the interac- 
tion one usually uses a different parameter kpa, where 



fc^/(6^^ 



Nf/V, and thus the relation between the 



two Fermi momenta is fc/ = kp{l 



)l/3 



III. Induced P-wave pairing gaps in homogeneous 
systems 

In Fig. [T]we show the P-wave pairing gaps computed 
using various approximations, namely the pairing gap 
computed in the BCS approximation used in Refs. (y, |9|, 
the pairing gap computed in the simple Eliashberg ap- 
proximation and finally the pairing gap computed 



using the full momentum and energy dependence as de- 
scribed above. In the case of P-wave pairing gap one 
has to distinguish between two possible angular quan- 
tum numbers, namely {l,m) — (1,0) when Aio(p) oc Pz 
and the pairing gap has a nodal plane, and (Z, m) = (1, 1) 
when Aii(p) px + ipy and the time reversal invariance 
is broken. In the BCS approximation the pairing gap 
has the same value in both cases. That is not the case ei- 
ther in the simple Eliashberg approximation or when the 
full energy-momentum dependence of the interaction is 
taken into account, and it appears that when full energy- 
momentum dependence is taken into account An > Aio, 
unlike the case of the simple Eliashberg approximation. 

The BCS approximation proves to be a rather inaccu- 
rate theory. The simple Eliashberg approximation under- 
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estimates the magnitude of the pairing gap by a signif- 
icant factor, sometimes close to an order of magnitude. 
In hindsight all this does not come as a surprise, since 
for these systems, both magnitudes of the exchanged en- 
ergy and momentum between two interacting fermions 
are large, and thus both effects due to a large range of 
the interaction and to retardation are naturally large. It 
is also natural to find that the magnitude of the pairing 
gap is a monotonous function of the relative dimer den- 
sity y — Hh/nf. The dependence of the pairing gap on the 
scattering length a is opposite to the dependence of the 
S-wave pairing gap in the symmetric phase. On the BEG 
side of the Feshbach resonance the pairing gap increases 
towards the BEG limit and tends to half of the binding 
energy of the dimer. The P-wave gap however shows the 
opposite behavior, it decreases in magnitude going from 
the Feshbach resonance towards the BEG limit (in which 
case the scattering length a decreases). The fact that 
the S-wave and the P-wave pairing gap have an opposite 
behavior is noteworthy. Another very important feature, 
which has already been noticed in Ref. is that the 
magnitude of the pairing gap becomes comparable with 
the Fermi energy of the surplus fermions when approach- 
ing the Feshbach resonance. Our calculation does not 
include the vertex corrections however, which, as many 
authors have argued are not expected to result in 
major changes. 

IV. Aspects of induced P-wave pairing in atomic 
traps and final remarks 

Since the induced P-wave pairing gap can be paramet- 
rically large while approaching the Feshbach resonance, 
the prospects to put it in evidence experimentally appear 
good. However, in current realizations of such systems in 
atomic traps, the role of the confining potential is going 
to be significant, and the realization of a homogeneous 
phase is unlikely. In a 3-dimensional trap one expects in 
the local density approximation a shell structure for an 
asymmetric system. The center of the trap will be occu- 
pied by a pure paired unpolarized phase, followed by a 
mantle where both dimers and surplus atoms will coex- 
ist, and at the outskirts a pure unpolarized phase [20| . 
It was shown in Ref. that in the weak coupling limit 
(when a — > -|-0) and in a harmonic trapping potential the 
radii of these shells are approximately related 

6Rl^Ri + 5Rl (8) 

where i?i is the radius of the pure dimer/unpolarized 
phase (x(r < i?i) = 1 and nf{r < = 0), i?2 is the 
radius where the dimer number density vanishes {nb(r > 
R2) = and x{r > R2) = 0), and i?3 is the radius of 
the cloud {nf{R^) = 0), thus Ri < R2 < R3. Dimer and 
surplus majority atoms coexist within the shell -Ri < r < 
R2. One can expect the appearance of a P-wave pairing 



only in this shell Ri < r < R2 where the ratio y = ni,/nf 
changes from y{Ri) = cxd to y(i?2) = 0, while e/(i?i) = 
and has a finite value at i?2. Thus, assuming that the 
local density approximation is valid, we see that there 
are two competing factors in determining the magnitude 
of the P-wave pairing gap. The local Fermi energy £/(?'), 
which determines the order of magnitude of the P-wave 
pairing gap, increase from zero at i?i to a final value at 
R2 

ef{r) (x5r^ + RI-6RI, (9) 

while the dimer to fermion ratio has an opposite behav- 
ior, decreasing from infinity to zero 

and consequently the P-wave pairing gap will attain a 
maximum absolute value in the middle of this shell. It 
will be a significant challenge for experimentalists to cre- 
ate a trap with a spatially rather large shell Ri < r < R2 
and a significant amount of surplus fermions so as to 
make visible this new exotic and quite unique in its na- 
ture P-wave paired phase. It is important to remember 
that kfa — kpa{l — x)^^^ = kpa/ {X + y^/'^ < kpa, where 
in a trap kj and kp are the local values. When one char- 
acterizes the proximity to the Feshbach resonance of a 
given Fermi gas, one typically uses the local Fermi mo- 
mentum for the majority component at the center of the 
trap, which is significantly larger than the local kp (and 
kf) in the region i?i < r < i?2 of interest here. We 
expect that for finite values of the scattering length the 
radii i?i,2,3 will change, but these changes will not be 
qualitative and the conclusions drawn by us above about 
the structure of the mixed atom-dimer shell Ri < r < R2 
remain qualitatively correct. 

In conclusion, we have determined the induced P-wave 
paring gap in an unbalanced system of atoms, where 
dimers coexist with atoms. Since the interaction be- 
tween the like fermions is mediated by sound modes and 
it has a strong energy and momentum dependence, it 
could be either attractive or repulsive in character, de- 
pending on the particular values of the exchanged energy 
and momentum between the two atoms. We have com- 
pared the BCS induced P-wave pairing gap with pairing 
gap determined in the simple Eliashberg approximation, 
when only the energy dependence of the interaction is 
taken into account, and in the full Eliashberg approxi- 
mation, when both momentum and energy dependence 
of the induced interaction is accounted for. We have ob- 
served that both the BCS approximation and the simple 
Eliashberg approximation, with only a frequency depen- 
dence taken into account, do not lead to correct values 
of the pairing gap. Moreover, we have observed that 
by approaching the Feshbach resonance the pairing gap 
increases substantially and thus the conditions that in- 
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duced P-wave pairing could be observed in atomic trap 
experiments improve significantly. 
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